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ABSTRACT 

Gamma-ray bursts (GRBs) are the most luminous astrophysical events observed so far. They 
are conventionally classified into long and short ones depending on their time duration, Tgo. 
Because of the advantage their high rcdshifts offer, many efforts have been made to apply GRBs 
to cosmology. The key to this is to find correlations between some measurable properties of 
GRBs and the energy or the luminosity of GRBs. These correlations are usually referred to as 
luminosity relations and are helpful in understanding the GRBs themselves. In this paper, we 
explored such correlations in the X-ray emission of GRBs. The X-ray emission of GRBs observed 
by Swift has the exponential functional form in the prompt phase and relaxes to a power-law 
decay at time T p . We have assumed a linear relation between logLx,p (with Lx, P being the 
X-ray luminosity at T p ) and log[T p /(l + z)], but there is some evidence for curvature in the data 
and the true relationship between Lx, P and T p j (1 + z) may be a broken power law. The limited 
GRB sample used in our analysis is still not sufficient for us to conclude whether the break is 
real or just an illusion caused by outliers. We considered both cases in our analysis and discussed 
the implications of the luminosity relation, especially on the time duration of GRBs and their 
classification. 



Subject headings: gamma-ray burst: general 
1. Introduction 

Gamma-ray bursts (GRBs), which can last 
from milliseconds to nearly an hour, are the most 
luminous astrophysical events observed so far. 
The parameter of T go , which is defined as the time 
interval during which the background-subtracted 
cumulative counts increase from 5% to 95%, is 
usually used to denote the time duration of GRBs. 
Those with T 90 > 2 s are conventionally described 
as long/soft GRBs and those with Tgo < 2 s as 
short /hard GRBs (jKouveliotou et allll993h . 

In lWillingale et all (|2007t ) . it was demonstrated 
that the X-ray decay curves of GRBs measured by 
Swift can be fitted using one or two components — 
the prompt component and the optional after- 



glow component — both of which have exactly 
the same functional form comprised of an early 
falling exponential phase and a following power- 
law decay. The prompt component contains the 
prompt gamma-ray emission and the initial X- 
ray decay. The transition time T p from the ex- 
ponential phase to the power-law decay in the 
prompt component defines an alternative esti- 
mate of the GRB durati on, which is compara ble 
with T 9n for most GRBs (I O'Brien et "aj] |2006aHbh 



O'Brien fc Willingald (|2007l ) proposed the clas- 



sification of GRBs into long and short ones at 



T p = 10 s instead. 



GRBs can be observed at very high redshifts 
due to their high luminosities. For example, the 
recently observed GRB 090423 has a redshift 
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of z » 8.2 (|Tanvir et all 120091: ISalvaterra et al 
2009). It may be possible to calibrate GRBs 



Lamb et al 



2005 



2005 



as standard candles (see, for example. iDai et al 
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Qi et al.ll2 008a.b; Kodama ct al. 2008; L iang et al 
Wandl2008l:lQi et alj|2009l etc.). The key to 



the calibration of GRBs is to find correlations be- 
tween some measurable properties (the luminosity 
indicators) of GRBs and the energy (the isotropic 
energy Ej^ so or the collimation-corrected energy 
Ej) or the luminosity (e.g., the peak luminosity 
L) of GRBs. These correlations are usually re- 
ferred to as luminosity relations, which are useful 
both in applying GRBs to cosmology and under- 
standing GRBs themselves. The GRB luminosity 
relations used in cosmological studies in the litera- 
ture include the relations of ri ag (the spectral lag, 
i.e., the ti me shift between th e hard and soft light 
curves)-L ( Norris et al. 2000h . V (the variability, 
a quantitative measurement on the spikiness of 
the light curve; there exist several definitions of V, 
mainly depending on the smoothing time intervals 
the reference curve is built upon, and on the nor- 
malization as well)-L (Fenimore fc Ramirez- Ruid 



20001 : iReichart et al.ll200lh. (peak e nergy of 

the spectrurcQ-E^isn dAmati et al. 2002), E VP ^ 



E^ (IGhirlanda et al 



2004al) . E peak -L (jSchaefer 



2003), an d trt (minimu m rise time of the light 



curve)-!/ (pchaeferl 120070 . For each of the above 
luminosity relations, there is only one luminosity 
indicator involved. More complicated luminosity 
relations, which include two luminosity indicators 
are also discussed in the literature; see, for exam- 
ple, |Yu_et_al. I (I2009Q and references therein. 

In this paper, we explore the correlation be- 
tween T p and the X-ray luminosity of GRBs at T p 
and discuss its implications, especially on the time 
duration of GRBs and their classification. 

2. Methodology 



In IWillingale et al 



(|2007l ). the X-ray light 
curves of GRBs were constructed from the com- 
bination of Burst Alert Telescope (BAT) and X- 
ra y telescope (XRT) da ta in the way described 
bv lO'Brien et all (|2006bl . the BAT data is extrap- 



olated to the XRT band) and fitted using one 
or two components — the prompt component and 
the optional afterglow component — both of which 
have the same functional form: 



fc(t) 



F c exp (a c - ^) exp , t < T e , 



t > T r , 



(1) 

The transition from the exponential to the power 
law occurs at the point (T c , F c ). The subscript c is 
replaced by p for the prompt component and by a 
for the afterglow component, and the fitted X-ray 
flux is the sum of the two components, i.e., f(t) = 
f p (t) + f a (t). In the derivation of the parameters, 
an initial fit was done to find the peak position of 
the prompt emission; the peak time was then used 
as time zero. A second fit was done with t p = 0. 
Thus, the derived parameters are all referenced 
with respect to the estimated peak time rather 
than the somewhat arbitrary BAT trigger time. 
In addition, large flares ha ve been masked ou t of 
the fitting procedure (see Lazzati et al. 2008L for 
a discussion on the average luminosity of X-ray 
flares as a function of time). 

We investigated the correlation between the 
transition time T p for the prompt component and 
the X-ray luminosi ty of GRBs at T p (see also 
Dainotti et al. 20081 for a discussion on the corre- 
lation between T a and the luminosity at T a ). The 
fit of the GRB light curves directly gives the values 
of T p . The X-ray luminosity of GRBs at a given 
time t is calculated by 



L x {t) = 4nd 2 L (z)F x (t), 



(2) 



where <Il{z) is the luminosity distance, for which 
we have used the flat ACDM cosmological model 
with fl m = 0.27 and H = 71 kms" 1 Mpc~\ and 
Fx (t) is the if -corrected flux given by 



F x (t) = f(t)x 
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f(t)x(l + zf-\ 



(3) 



where j3 is the spec tral index and, in gene ral, it 
is time dependent. IWillingale et al.l (|2007| ) pre- 
sented in Table 3 of their paper the spectral index 
in the prompt phase (/3 P ; from the BAT data), in 
the prompt decay (f3 p d', from the XRT data), on 
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the plateau of the afterglow component (/3 a ; XRT 
data), and in the final decay (Pad', XRT data). 
Since the BAT data have been extrapolated to the 
XRT band in the combination of the BAT data 
and XRT data, {E min , £ max ) = (0.3, 10) keV 
should be used here and for this limited energy 
range, the simple power-law is sufficient for 
the GRB spectrum. So, the X-ray luminosity of 
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GRBs at the transition time T p is 



And 



l( z ) f x, p 
(z)F p (l + zf"- 



(4) 



where we have used (3 — f3 p at t = T p . 

For the investigation of the correlation, we 
mainly fit the data to the relation 



log L x „ 



b log 



T„ 



1 



(5) 



where T p /(1 + z) is the corresponding transition 
time in the burst frame. Like many other lumi- 
nosity relations, this relation is by no means an 
accurate one. Due to the complexity of GRBs, 
it is hardly possible that \ogLx, P is fully deter- 
mined by only log[T p /(l + z)}. As usual, intrin- 
sic scatter <7i nt is introduced here, i.e., lacking 
further knowledge, an extra variability that fol- 
lows the normal distribution A/"(0, af nt ) is added 
to take into account the contributions to log Lx, P 
from hidden variables. For the fit, we used the 



techniques presented in D'Agostini ( 20051) . Let 
x = log[T„/(l + z)} and y = \ogLx, P ; according 



to lD'Agostinil (|2005j ). the joint likelihood function 
for the coefficients a and b and the intrinsic scatter 
o-int is 



C(a, b, (Tint) oc Y[ 



x exp 



+ < + b2 < 
(y l - a - bXiY 



2 R t 



6 V? 



(6) 



where xi and yi are corresponding observational 
data for the zth GRB. 

For G RB data, we u se d the samples com- 
piled in IWillingale et all (|2007t) . that is, the 
107 GRBs detected by both BAT and XRT 
on Swift up to 2006 August 1. The T 90 du- 
ration for these 107 GRBs were obtained from 



ties o f the fitted parameters in IWillingale et al 
(|2007l ) are given in 90% confidence level. We sym- 
metrize the errors and derive the corresponding la 
uncertainties by just dividing the 90% confidence 
level errors by 1.645. Unless stated explicitly, 
the errors in this paper are for the la confidence 
level. In our analysis, when not all of the GRBs 
have the needed parameters available, we use the 
maximum subset of the 107 GRBs satisfying the 
requirement. For example, in Eq. (j4]), the calcu- 
lation of Lx,p needs the observed redshift z, in 
addition to F p and (3 P , which are derived from the 
fit of GRB data. 

3. Results and discussion 

We plot the logarithm of Lx, P versus the log- 
arithm of T p /(1 + z) in Figure [TJ which include 
47 GRBs. We can see that most of the GRBs 
(34 GRBs) lie in the range of 2 s < T p /(1 + z) < 
100s. There is obviously a correlation between 
Lx,p and T p /(1 + z). However, when fitting the 
GRBs to the relations of Eq. ([5]), we have differ- 
ent options depending on how we view the three 
GRBs with the largest T p /(1 + z), i.e., those with 
T p /(1 + z) > 100 s. For the first choice, we can 
simply include all the data points in the fit (see 
the top panel of Figure [2]) , which leads to a re- 
sult that the GRB with the largest T p /(1 + z) lies 
outside the 2a confidence region of the fit. Alter- 
natively, it is also possible that the GRBs with the 
largest T p /(1 + z), instead of just being outliers, 
may indeed reveal some trend of the luminosity 
relation at large T p /(1 + z). In this case, we can- 
not simply ignore the GRBs with large T p /(1 + z) 
just because their quantity is small and, if they 
are taken more seriously, it seems that the sam- 
ples are split into two groups at some value of 
T p /(1 + z) based on the slope of the luminosity 
relation of Eq. ([5]). To show this, we perform a 
fit using only the GRBs with T p j (1 + z) > 2 s (see 
the bottom panel of Figure [5]), which gives a result 
quite different from the prior fit. From the bottom 
panel of Figure [21 we can see that if the best fit 
line is extended to the range oiT p /(l + z) < 2s, 
all the GRBs with T p /(1 + z) < 2 s lie below the 
line. As a comparison, we also fit the GRBs with 
2s < T p /{l + z) < 100s (see Figure[3]), which show 
that the difference is indeed introduced by the 
three GRBs with the largest T p /(l + z) when com- 



3 



pared with the fit of GRBs with T p /(1 + z) > 2 s. 
We tabulate the fit results in Table [TJ From the 
table, we can see that the values of b for the first 
two cases are quite different. In addition, it is 
also interesting to note that, for the first case, 
the slope b is close to the slope (— O^l^g) of a 
similar luminosity rela tion about T a and lum inos- 
ity at T a presented in iDainotti et al.l ( 20081 ). For 
the second case, the slope b is close to the index 
(—1.5 ± 0.16) of the power-law declination of the 
average luminosity of X-ray flar e s as a function of 



(|2008fi . These co- 



time presented in lLazzati et al. 
incidences may be worth noting in future studies 
with a bigger sample of GRBs. 



GRB sot 










b 




"int 




All available GRBs 


50 


91 ± 


23 


-0 


89 ± 0.19 


1 


06 ± 


13 


T P > 2s 


51 


96 ± 


32 


-1 


73 ± 0.25 





78 ± 


11 


2s < < 100s 


51 


09 ± 


32 


-0 


74 ± 0.30 





63 ± 


09 




Table 1: Results of the fit to the luminosity rela- 
tion of Eq. ©. 




Fig. 1. — Lx, P (in ergs -1 ) versus T p /(1 + z) (in 
seconds). 47 GRBs are included with 37 GRBs 
that have T p /(1 + z) greater than 2 seconds and 3 
of which have T p j (1 + z) greater than 100 seconds. 
The red ones are conventional short GRBs (Tg < 
2s). 

Generally speaking, in statistical analysis, ex- 
cept for the different measurement precisions, we 
should treat all the data points equally and should 
not give more attention to some data points over 
the others. However, since our sample of GRBs 
here is very limited, and there are only a few GRBs 
with very large T p /(1 + z), some selection rules 



Fig. 2. — Best fit of the GRBs to the luminosity 
relation of Eq. ([3]) and the corresponding 2a con- 
fidence region. Top: all GRBs are included in the 
fit. Bottom: only those with T p /(1 + z) > 2 s are 
included in the fit. 
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Fig. 3. — Best fit of the GRBs to the luminosity 
relation of Eq. ([5|) and the corresponding 2a con- 
fidence region. Same as Figure [3] except that only 
GRBs with 2 s < T p /(1 + z) < 100 s are included 
in the fit. 



may have already been imposed implicitly on the 
sample itself. Considering this, it may be unfair 
to the GRBs with large T p /(1 + z) to be treated 
as outliers just because their quantity is small, es- 
pecially when they may reveal some trend in the 
luminosity relation. This is why we consider both 
cases in the analysis, i.e., whether the few GRBs 
with large T p /(1 + z) should be treated as just 
outliers or taken more seriously. Correspondingly, 
the relation between Lx, P and T p j (1 + z) could be 
a simple power law or a broken power law with a 
change in the slope of Eq. §5§ at some character- 
istic value of T p /(1 + z). For the present sample 
of GRBs used in our analysis, it is not sufficient 
for us to conclude which case is real and, for the 
later one, to determine the exact value of T p J (l+z) 
where the slope of the luminosity relation changes. 
Here, we leave it open to future studies with more 
GRBs and discuss in the following the implications 
of the luminosity relation, especially in the situa- 
tion where there is a change in the slope at some 
value of T p /(l + z). 

First of all, we emphasize that the luminosity 
relation is between the luminosity and T pi though 
T p and Tgo can both act as an estimate of the 
GRB duration and are comparable to each other 
for most GRBs, as can be seen from Figure 0] A 
similar relation seems not to exist between the lu- 
minosity and the Tgo. See Figure [SJ the corre- 



sponding data points turn out to be very disper- 
sive. Despite the fact that most GRBs in Figure |4] 
are distributed around the line on which they are 
equal, there are some of them that have been con- 
siderably different from T p and Tgo. In fact, the 
two quantities differ from each other significantly 
from their derivation. Tgo is calculated directly 
by using the BAT data in the 15-150 keV band, 
while for the calculation of T p , the BAT data are 
first extrapolated to the XRT band of 0.3-10 keV 
in order to be combined with the XRT data. For 
Tgo, the emphasis is on the percent of the fluence 
of a burst, while for T p , the emphasis is on the 
transition in a GRB light curve from the expo- 
nential decay in the prompt phase to the initial 
power- law decay. In the time interval of Tgo, 90% 
of the total fluence is observed, while the ratio be- 
tween the observed fluence in the time interval of 
T p and the total fluence depends not only on the 
temporal decay index of the initial power-law de- 
cay duration in the prompt component, but also 
on the shape of the light curve in the afterglow 
component. In addition, we must remember that 
large flares have been masked out from the light 
curves bef ore performing the fit that allows us to 
derive T p (jWillingale et al.ll2007l ). 




Fig. 4.— T p versus T 90 for the 107 GRBs. Times 
are in seconds. On the cyan dash-dotted line 
the two quantities have equal values and the red 
dashed lines correspond to T p and Tgo equal to 2 
seconds respectively. 

As stated above, because the data are limited, 
there are two possible models for the luminosity 
relation of Eq. ([5]), i.e., the luminosity relation, 
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Fig. 5. — Lx, v (in ergs^ 1 ) versus T 90 /(l + z) (in 
seconds). The set of GRBs is the same as that 
in Figure [TJ The red ones are conventional short 
GRBs (T 90 < 2 s). 

using one set of values for the parameters (the in- 
tercept a and the slope b), may be applicable to 
all the GRBs except for some outliers, or there 
may be a change in the slope b of the luminosity 
relation at some value of T p /(1 + z). If there is a 
change in slope this may suggest that GRBs could 
be classified into two groups based on their val- 
ues of T p / (1 + z). Since T v j (1 + z) is an estimate 
of the GRB duration, this is in fact an indica- 
tion of how we should classify GRBs into long and 
short on es and is actually the same as the pro- 
posal bv lO'Brien fc Willingalel (|2007t ) to use T p as 
the criterion for the classification of long and short 
GRBs. In principle, we should use the quantity in 
the burst frame (T p /(1 + z)) instead of that in the 
observer frame (T p ). However, for a large portion 
of the observed GRBs, the redshifts are not avail- 
able. Generally speaking, due to the diversity of 
GRBs, the classification of a GRB is unlikely to be 
completely determined by only its time duration, 
not to mention that the time duration of long and 
short GRBs overlaps near the demarcation point. 
Let us assume the demarcation point in the burst 
frame for the long and short GRBs to be 



T p /(1 + z) = T l/S 



(7) 



Then, as an approximate method, an effective red- 
shift z c ff can be defined, such that GRBs can be 
classified into long and short ones in the observer 
frame at 



In addition, since a similar relation does not hold 
if we replace T p with Tgo as mentioned previously, 
the change in the slope b, which may be a reflection 
of different mechanisms, if confirmed, would favor 
T p over Tgo as a criterion for the classification of 
long and short GRBs. 

4. Summary 

We investigated the correlation between T p and 
the X-ray luminosity of GRBs at T p and found 
a (broken) linear relation between \ogLx, P and 
log[T p /(l + z)}. There may be a change in the 
slope of the relation at some value of T p /(1 + z) 
mainly because of the presence of the few GRBs 
with large T p /(\ + z). The limited GRB sample 
used in our analysis is still not sufficient for us to 
conclude whether the change in the slope is real 
or just an illusion caused by outliers. We consid- 
ered both the cases in our analysis. If the change 
is real, the different slopes may be a reflection of 
different mechanisms for GRBs, which may sug- 
gest that using T p instead of Tgo (considering that 
a similar relation does not hold if we replace T p 
with Tgo, though Tgo and T p are both estimates of 
the GRB duration) as a criterion for the classifi- 
cation of long and short GRBs. 
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